arXiv:1505.01561v2 [cond-mat.mes-hall] 26 Oct 2015 


Up-converted photoluminescence induced by radiative coupling between excitons 

Takuya Matsuda0 Nobuhiko Yokoshi, and Hajime Ishihars0 
Department of Physics and Electronics, Osaka Prefecture University, Sakai, Osaka 599-8531, Japan 

We propose an unconventional scheme of photoluminescence in a semiconductor thin Him, where 
the nonlocal correlation between an excitonic wave and a light wave prominently enhances the 
exciton -light coupling beyond the long-wavelength approximation (the so-called excitonic superra¬ 
diance regime). On the basis of the developed method extending input-output theory, we elucidate 
atypical photoluminescence effects due to the strong wave-wave correlation. In particular, the up- 
converted photoluminescence based on the coherent superposition of excitons is found to be highly 
efficient, i.e., it can be realized by weak pumping without auxiliary systems such as cavities or 
photonic antennas. 

PACS numbers: 78.67.-n, 71.36.+C, 42.50.Nn 


Photoluminescence (PL) is one of the fundamental 
phenomena arising from light-matter coupling, and has 
been extensively utilized for monitoring the structure 
of the electronic levels in materials 1]. PL spec¬ 
troscopy provides the essential information of matter sys¬ 
tems, such as the structure of the levels and oscillator 
strengths of electronic excited states after reaching quasi- 
equilibrium. On the other hand, PL is also important for 
observing the interplay scheme between light and elec¬ 
tronic systems. A representative example is PL reflecting 
the dispersion relation of polaritons and not the thermal 
equilibrium for bare excitons SHU- In order to access 
polaritons by photon emission, photons have to be emit¬ 
ted before the exciton-photon coupling is “dephased”, 
or the excitations relax to quasi-equilibrium. Then, it is 
necessary to consider polariton accumulation at the bot¬ 
tleneck of the decay SI, or confine the electronic system 
in high-density photons @-@|. 

Recently, another scheme that allows PL to access 
exciton-light coupled modes has been demonstrated in 
simple semiconductor thin films @. When the film thick¬ 
ness is in the nano-to-bulk crossover regime (a few hun¬ 
dred nanometers), the exciton light coupling is extremely 
enhanced by the nonlocal correlation between a light 
wave and an excitonic center-of-mass (c.m.) wave [lOj— 
fl3j |. Note that the resultant exciton light coupled state 
should be distinguished from polaritons because it ex¬ 
hibits an ultrafast radiative decay exceeding the polari¬ 
ton formation, and is in the so-called excitonic superra¬ 
diance regime fT3 - fl6| . Such a strong coupling beyond 
the long-wavelength approximation (LWA) leads to un¬ 
precedented PL properties. Actually, photon emission 
simultaneously occurs from more than one exciton light 
coupled mode including optically-forbidden ones @. 

Here, we focus on a unique feature of the beyond the 
LWA regime; there exists strong coupling between the 
different c.m. states of excitons via radiation that causes 
quantum superposition of these states. The purpose of 
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this paper is to theoretically demonstrate that the control 
of the superposition provides unconventional and func¬ 
tional aspects of PL properties. In particular, remarkable 
up-converted PL is revealed, where the up-conversion ef¬ 
ficiency is greatly strengthened for appropriate system 
parameters. Up-converted PL has been an appealing 
subject not only for its fundamental interest but also 
for i ts p otential applications such as frequency conver¬ 
sion 1171123l |. There are some typical schemes, e.g., those 
based on phonon-assisted processes [l7l - [20| . and on mul¬ 
tiphoton processes [21N23I ] . Recently, Fernee et al. re¬ 
ported that the up-conversion range is greatly widened 
by an unusually-large dephasing rate in a colloidal so¬ 
lution of quantum dots [241 . The up-conversion in the 
present work is essentially distinguished from the pre¬ 
vious ones in that it is the PL from the exciton- light 
coupled state with a higher exciton light coupled mode. 
The appearance of such a higher state can be explained 
as follows; a portion of the excited eigenmode is recon¬ 
structed to superposed eigenmodes with dephasing into 
different c.m. exciton states. In addition, the main pro¬ 
cess can be clearly realized by weak pumping without 
auxiliary systems such as cavities or photonic antennas. 

In the following demonstrations, we treat an explicit 
geometric model of thin films for numerical calculations 
in order to discuss the observed results in Ref. H The 
thin film confines the c.m. motion of an exciton per¬ 
pendicular to the surface and acts as a homogeneous, 
nondispersive, and nonabsorptive background medium 
with a dielectric constant £b g for the excitons. For de¬ 
veloping the theoretical method to calculate the lumines¬ 
cence spectra for the nano-to-bulk crossover size regime, 
it is necessary to take account of the nonlocal optical re¬ 
sponse of the excitons; thus, we explicitly consider the 
microscopic spatial structures of the excitonic wave and 
light wave. The Bohr radius of an exciton is assumed to 
be much smaller than the film thickness L; therefore, its 
relative wavefunction is approximated to be the same as 
that in a bulk system. Thus, neglecting the distortions 
near the surfaces, we assume the c.m. states of exciton 
to be simple sinusoidal waves: = 2/Zv sin(fc M z), 

where k^ = [in/L is the quantized wavevector with 
[i = 1, 2, • • •. The transverse energies and the eigenen- 
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ergies of the exciton are given as fiwx = 3.2022 eV and 
hfi e * = Huj t + (Hk^) 2 /(2m ex ), respectively. Here, m ex is 
the translational mass of the exciton. 

The Hamiltonian of the exciton- light coupled system 
is written as H = H ex + H Tai j + H lnt . The Hamilto¬ 
nian of the excitons is described as H ex = M , 

and b,j, (fef) stand for the bosonic annihilation (creation) 
operator of the /xth exciton state. The Hamiltonian of 
the radiation field is represented as IT ra d = Ej? 
where a v (al) means the annihilation (creation) opera¬ 
tor of the 77 th photon mode with an energy hfi v . The 
interaction between the exciton and the radiation field 
is expressed as H lnt = — f dzP ex (z)E(z). Here, the 
excitonic polarization operator P ex (z) is represented as 
Ax(z) = E M ( / P^{z)b li +U.c.) by using V^z) = V ip^z) 
with the dipole moment V. The moment V can be 
estimated by the longitudinal-transverse (LT) splitting 
energy Tzwlt = T’ 2 /(eo£bg)- The operator E{z) = 
a v (i£T)(z)ari + H.c.) denotes the electric field, where 

a v = yjhfirj /( 2 ero), £0 is the vacuum permittivity and 
£ v (z) is the eigenfunction satisfying Maxwell equation 
without excitonic polarization sources. 

As for the excitonic non-radiative processes, we treat 
them under the Born-Markov approximation. Then, the 
density matrix of the entire system obeys the following 
master equation: 


d 1 - 

-Qtm) + Aamp/5(f) + ^phaseP(i), (1) 


where the non-radiative decay of the exciton states, 
e.g., due to phonons, is described by £damp/S(i) = 

E M (7ex/ 2 )[ 2 0'M.M+lP( i )^ ,/i+l P(t)} 

and the dephasing is described by £ p ha.seP{t) = 
Eu(- r ex/2) fyfip, [hlpp, p(t)] . We have introduced the 


lowering operator for excitons cr M . M+ i = b^b^+i. In the 
present demonstration, it is assumed that non-radiative 
transitions only occur between the nearest levels. All 
non-radiative decay rates are the same tvy ex = 0.01 meV, 
and the dephasing rates of each exciton state are hr ex = 
0.2 meV considering the cryogenic conditions. These as¬ 
sumptions do not affect the essence of the following dis¬ 
cussion because the non-radiative decay rates are much 
smaller than the radiative decay rates under the present 
conditions, although the quantitative details of the re¬ 
sults depend on the values of the non-radiative ones (the 
influence of Ay ex is examined in Ref. [Ell . From Eq. ©, 
we obtain the simultaneous Heisenberg equations to de¬ 


termine the polarization ( 6 ^), the excitonic population 
( b^bp), and the correlation f25|. The detailed so¬ 

lutions of these equations are given in Ref. [26|. 

First, we analyze the excitonic population (b^bf,). 
Hereafter, unless otherwise noted, the system is assumed 
to be a CuCl film with a thickness L = 325 nm model¬ 
ing the experiment in Ref. J^j, where the dielectric con¬ 
stant is £b g = 5.59. As for the incident light, we em- 
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FIG. 1: (Color online) (a) The dotted lines represent the 
dispersion curves of the upper- and lower-branch polaritons 
(UBP and LBP), and the cross marks indicate the resonance 
energies of the exciton light coupled modes (vertical dotted 
lines) in the film with L = 325 nm. The bar length represents 
the radiative decay rates of the coupled modes, m denotes the 
index of the original exciton state of each coupled mode. (b,c) 
The population spectra of different exciton states (in this case 
= 1, 2, • • • ,8) are plotted versus the incident energy, and fj, 
represents the exciton state. Every population spectrum has 
a complicated structure because of the very large radiative 
coupling between different excitonic waves. All populations 
are normalized by the peak value of the population of the ex¬ 
citon states (til. The c.m. wavefunctions of the excitons 
[see inset of (b)] are considered. 


ploy a continuous wave laser whose intensity is A n = 
C£ P £ 2 n /2 = 100 W/cm 2 , where c is the speed of light, and 
fin is the amplitude of the incident light. In Fig. HJa), 
we plot the eigenenergies of the exciton light coupled 
modes, as denoted by the cross marks, whereas the dot¬ 
ted lines indicate the dispersion relation of the upper- 
and lower-branch polaritons. It should be noted that the 
dispersion relation of the exciton- light coupled modes 
(to = 1 , 2 ,...) deviates from that of the polaritons in 
the bulk system owing to the finite radiative width a 
The radiative width becomes enlarged near the LT split¬ 
tings, and it reaches 96 meV for the coupled mode to = 5. 
Figures [U[b,c) show the populations of the different ex¬ 
citonic states plotted versus the incident light en¬ 

ergy. One can see that each population has peaks reflect¬ 
ing the eigenenergies of the exciton light coupled modes 
(vertical dotted lines) due to the radiation-mediated cou¬ 
pling between the different c.m. states of the exciton. 
This coupling is described with the terms appearing in 
the Heisenberg equations (see, El): 

Zw> =~J dz J dz ' WGMMA ( 2 ) 

where G(z,z') is the Green’s function for Maxwell equa- 
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tion in the vacuum/film/vacuum structure j27j ]. This 
term represents the coupling between the exciton states 
mediated by the radiation fields. Because of this inter¬ 
action, each population (b^b^) exhibits structures due 
to the resonances at the exciton light coupled modes or 
the interference between them. Notice that the peaks 
slightly deviate from the eigenenergies especially above 
the LT splitting, because the adjacent modes overlap ow¬ 
ing to the radiative width. By focusing on the c.m. states 
fi = 4, 5, we can clearly see that their spectra remain vis¬ 
ible on both the sides of the LT splitting. This means 
that the wave-wave coupling between the excitons and 
the radiation field causes the coherent superposition of 
different exciton c.m. states over the LT splitting. 

Subsequently, we introduce a calculation method to in¬ 
vestigate the PL properties. Here we extend the conven¬ 
tional input-output theory to include the nonlocal cor¬ 
relation inside the film. By using the Green’s function 
for the radiation field, the electric field can be 
written as (28j 


E(z 0 ,t) = E 0 {z 0 ,t) + y; f dz' G(z 0 ,z')V IJ .{z')b^{t), 

(3) 


where zp is an arbitrary position. Here the incident field 
Ep{zp,t) corresponds to the one in the absence of the 
resonant contributions of the excitons. Then, the in¬ 
put pump field is the plane wave with the amplitude 
fin, which is connected to satisfy the Maxwell bound¬ 
ary condition in the vacuum/film/vacuum structure. By 
using the output radiation field in Eq. ©, we calculate 
the first-order correlation function (E\zp, t)E(zp, t + t)) 
according to the quantum regression theorem |29|, |30| . 
and decompose it into a coherent component and an 
incoherent one. The coherent component corresponds 
to elastic scattering, whereas the incoherent compo¬ 
nent provides inelastic scattering, namely “the lumi¬ 
nescence” . Then, the incoherent component is written 
as {AE\zp 1 t)AE{zp,t + r))={E^{zp,t)E{zp,t + r)) - 
{E^(zp,t)) (E(zo,t + t)), where the last term on the 
right-hand side corresponds to the coherent component. 
Under steady-state conditions, the PL spectrum on the 
transmission side is obtained by calculating the Fourier 
transform of the incoherent correlation function: 


SicM = Re 



e iuT - j. 

dr- (AE\,(zq, 0)AEt(zo, t)) 

7r 


(4) 


where the field ET(zp,t) shows the electric field on the 
transmission side. Hereafter we normalize the PL inten¬ 
sity by WLT/^in- 

In Fig. HJa), we plot the PL intensities versus the out¬ 
put energy for different incident light energies which are 
set to be resonant with the exciton light coupled modes 
below the LT splitting, i.e., m = {6,7,8} in Fig. Oja). 
One can see that luminescence beyond the LT splitting 
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FIG. 2: (a) The PL intensity is plotted versus the output pho¬ 
ton energy. The incident energy hu>i n is tuned to the resonance 
energies of the coupled modes m = 6, m = 7, and m = 8 in 
Fig.l (a), (b) The PL intensity is plotted under the condition 
in which the incident energy fku i n is tuned to 3.186 eV (exper¬ 
imental condition in Ref. |9j) and 3.203 eV (stopband region). 
The vertical dotted lines represent the eigenenergies of the 
exciton-light coupled modes. 



appears, i.e., the light energies are up-converted. Here, 
we remark that the energy conservation law should hold 
true for the whole system in steady state, even though the 
luminescence includes up-converted components. Indeed, 
the present model and calculation scheme guarantee this 
law, which is verified in Ref. [26j. 

The essential factor of the observable PL up-conversion 
is the exciton- light coupling, extremely enhanced by the 
nonlocal correlation beyond the LWA. The mechanism of 
PL up-conversion is interpreted to be a result of quantum 
superposition between a resonantly-excited exciton light 
coupled mode and the one above the LT splitting. To be 
precise, a small portion of the excited coupled mode is 
reconstructed into multiple superpositions of the coupled 
modes by the enhanced exciton-light coupling under slow 
dephasing to the exciton c.m. states (/i = 1,2,...). Ac¬ 
tually, the PL appears to reflect the eigenenergies of the 
coupled modes (vertical dotted lines), unlike conventional 
PL [lj. This means that the PL signal can plausibly iden¬ 
tify the exciton-light interactions. 

Further, photon emission is affected by the parity- 
selection rule of the exciton c.m. state. For example, 
if the coupled mode m =7 is resonantly excited, a signal 
predominantly appears from the odd-numbered states. 
Conversely, if the coupled mode to = {6,8} is excited, 
the even-numbered states provide the primary contribu¬ 
tion to PL. This is consistent with our criterion that 
the proposed up-converted PL originates from neither 
phonon assisted, multiphoton processes nor huge dephas¬ 
ing but from the radiation-mediated superposition of the 
exciton-light coupled modes. In addition, there must 
exist a finite probability amplitude, even inside the LT 
splitting, due to the overlap between the adjacent coupled 
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modes to superpose the states on both sides of the LT 
splitting. Indeed, the coherent superposition also affects 
the PL properties inside the LT splitting. In Fig [2jb) , 
the PL spectrum for hui i n = 3.203 eV (green dashed line) 
that lies inside the LT splitting. One can apparently see 
that PL signals are found inside the LT splitting, which 
never occurs in a bulk sample. 

In order to compare the calculation with the exist¬ 
ing experiment using a 389-nm coherent light source 
(Ti:sapphire) (T) . we investigate the case where the in¬ 
cident light energy is set to be huj in = 3.186 eV [see the 
red solid line in Fig. [U[b)]. Although the incident energy 
is far detuned from the coupled modes, we confirm that 
the signal is small but appears at the eigenenergies above 
the LT splitting. It should be remarked that we can actu¬ 
ally find a similar signal, namely, up-converted PL, in the 
experiment described in Ref. [9j [32|. Although the veri¬ 
fication of the origin of this signal by theory considering 
the detailed experimental conditions will be a subject of 
future study, this fact convinces us of the reality of the 
present proposal. As for the up-conversion efficiency, we 
evaluate the photon-number efficiency as 



Here, we count the photons emitted from the coupled 
state m = 4 above the LT splitting, i.e., between 
Huib = 3.220 eV and = 3.212 eV. Although the inci¬ 
dent intensity is as small as 100 W/cm 2 and off-resonant, 
the evaluated photon-number efficiency is vi' 2 i 2 eV = 
0.077%, which is sufficiently large to be detected in an 
experiment. 

Because the strong nonlocal correlation in the film re¬ 
quires synchronization between the excitonic wave and 
the light wave, the up-conversion efficiency should de¬ 
pend on the film thickness L. In Fig. [3|a), we show the 
PL intensity plotted versus the thickness and output en¬ 
ergy when 389-nm coherent laser light (huj\ n = 3.186eV) 
is injected. As is expected, the intensity is found to ex¬ 
hibit a strong thickness dependence. One may consider 
that the larger radiative width is favorable for efficient 
up-converted PL. However, when comparing Figs.[3][a,b), 
the peaks in the PL intensity do not directly coincide with 
those of the radiative width. We therefore investigate the 
excitation populations of the c.m. states of the exciton 
for the two thicknesses (see Fig. Etc)). One is the thick¬ 
ness used in Figs. |T| and [2 \{L = 325 nm), whereas the 
other is the thickness that provides a strong PL inten¬ 
sity (L = 280nm). By calculation, the photon-number 
efficiency for the thickness L = 280 nm is = 

0.207 %, which is more than twice as large as the one for 
L = 325 nm. By considering the reason for the enhance¬ 
ment, we now focus on the region around huji n = 3.19 eV, 
where the eigenenergy of the coupled mode m = 5 lies. 
One can see that the population of the exciton c.m. state 
fi = 5 has a broad peak for the thickness L = 280 nm (red 
solid line), whereas the very large radiative width satu¬ 
rates the peak for L = 325nm (pink dashed line). This 
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FIG. 3: (a) The PL intensity is plotted as functions of the 
film thickness and output photon energy. The incident energy 
is tuned to 3.186 eV. (b) The film thickness dependence of 
the calculated radiative width, (c) The population spectra of 
different exciton states (in this case fj, = 4, 5 at thicknesses 
of L = 280 nm, 325 nm) is plotted versus the incident energy 


leads to a larger population for the coupled mode m = 5 
at hu!i n = 3.186 eV for L = 280 nm, and an enhance¬ 
ment in up-converted PL. This means that efficient up- 
converted PL by off-resonant pumping can be obtained in 
the experiment by invoking the theoretical design of the 
system considering the populations of the exciton c.m. 
states. 

In summary, we have examined photoluminescence 
(PL) in a simple semiconductor thin film and proposed 
an unconventional type of up-converted PL in the exci¬ 
tonic superradiance regime. In order to treat the non¬ 
local correlation beyond the long-wavelength approxima¬ 
tion (LWA), we extend conventional input-output the¬ 
ory [29j, |30[. In our numerical calculations of the PL 
spectra, we have found that a portion of the incident en¬ 
ergy is up-converted, and the efficiency is enhanced by 
choosing the sample size and pumping frequency. The 
favorable conditions for enhancement clearly reflect the 
parity combination of the c.m. wavefunctions. The up- 
converted PL originates from the coherent superposition 
of different exciton- light coupled modes in the resonantly 
populated excitons. The strong exciton light coupling 
via the nonlocal correlation beyond the LWA greatly en¬ 
hances the superposed component energetically far from 
the excitation energy. This is why the main process can 
be clearly realized by weak pumping without auxiliary 
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systems such as cavities or photonic antennas. 

It should be noted that the thin film geometry used for 
the present demonstrations that exhibits a suggestive sig¬ 
nal in the experiment [§] is a very familiar system, and 
further, the proposed mechanism has generality. How¬ 
ever, the excitonic sup erradiance can be enhanced in pla¬ 
nar microcavities |33l l34|. Thus, it might be possible to 
design more efficient systems if one examines such sys¬ 
tems over a widely expanded volume [H, [3(|. By con¬ 


trolling the coherent light-matter coupling in designed 
structures, the proposed up-converted PL will lead to fu¬ 
ture studies of functional PL for novel light sources as 
well as sensitive probes of light-matter coupling. 

We thank M. Ashida and M. Ichimiya for valuable dis¬ 
cussions, and T. M. thanks H. Ajiki for fruitful discussion. 
This work was partially supported by a Grant-in-Aid for 
Scientific Research (A) No. 24244048 from Japan and by 
the Japan Society for Promotion of Science (JSPS). 
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Appendix A: DETAILED SOLUTIONS OF HEISENBERG EQUATIONS 

In this section, we present detailed solutions of the Heisenberg equations to obtain the expectation values of the 
polarization (6 M ), the population (bj^bfj) and the correlation (b^b^'} under steady-state conditions. In the main text, 
we consider the Hamiltonian as 


H = H ex + H iad + H int . (Al) 

The Hamiltonian H ex describes the excitonic system as 

He, = E (A2) 

where b M (b^) stand for the bosonic annihilation (creation) operator of an exciton state /i and 12® X is its eigenfrequency. 
The Hamiltonian iL ra( j represents the radiation field as 

H ra d = ^ CL^h'q, (A3) 

where a v (a 1) stand for the annihilation (creation) operator of the 77 th photon state with an energy hfi v . The coupling 
Hamiltonian H lnt signifies the interaction between the exciton and the radiation field as 

Hint = - J dz P ex (z)E(z). (A4) 

The excitonic polarization operator P ex {z) is represented as 

P ex (z) = + H.c.) (A5) 

by using the coefficient V^(z) = V ^^(z) with the dipole moment V. The electric field operator E(z) is denoted as 

E(z) = '^2a 71 {i£ r] {z)a ri +R.c.), (A6) 

where a v = y/hfi^/ (2eo) and £r](z) is the eigenfunction satisfying the Maxwell equation without excitonic polarization 
sources. By solving the Heisenberg equation of motion for the exciton and the radiation field, we can derive the electric 
field by using Green’s function for the radiation field [28| given by 

E(z,t) = E 0 (z,t) + J2 Uz’ G(z,z')V^z'%{t). (A7) 

We consider the density matrix of the whole system obeying the master equation in the main text given by 

d 1 r „ , , « 

= 7+ .Cdamp p(t) + .Cphas eP{t), 


(AS) 
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where the non-radiative decay of the exciton states is described as 


^dampp(t) — El 2 2CT AlljU +lp(i)<T^ !(U+1 {d'E+l^’MiAH-l ) P{t)} 
and the dephasing is given by 


4 >has eP(t)=J 2 ~T' KWW] 


2 L 


(A9) 


(A10) 


We introduce the lowering operator of the excitons dy^+i = bjE+i- We assume that the non-radiative transitions 
only occur between nearest levels. All non-radiative decay rates are the same 7 ex , and all dephasing rates are the 
same A ex . From the master equation in Eq. (IA8I) . we can obtain the following equations: 

(All) 


| <&„(*)> = -(i + 7-/2 + T ex /2) <6„(t)) + i f dz V;{z) (E{z, t )), 




- 7 ex(^(t) 6 M (t)) + 7 ex (t)V+l(t)) 

+ ^J dz {(t>l(t))V;(z)(E 0 (z,t))-V„(z) (El(z,t)) <S„(i)>) 

- { E ^ ^ E : 


(A12) 


"oT (bUt)b^(t)) = [i(12- - fl?) - 7ex - F ex l (bUt)b,,(t)) 


where 


dt 


+ - h J {(bl(t))r*,(z) (Eo{z,t))-V^(z){El(z,t)){b^{t))) 

- 4 E$+ 4E <&E)Eco>, 


= -J dz J dz ' K( Z ) G ( Z ’ Z '^'( Z ')- 


(A13) 


(A14) 


This term represents the radiation-mediated coupling between excitons, which plays a central role in our work. This 
radiative coupling between excitons causes the coherent superposition of different exciton states mediated by the 
strong exciton-photon couplin g vi a nonlocal correlation. By considering the self-consistent coupling between these 
equations and the electric field [2f|, the expectation values can be obtained. By substituting Eq. (1A71) , we can rewrite 
Eq. (IA11I) under the rotating wave approximation as 


E M m,a( w ) (VM) = J dz KW (Eo{z,u)) 


where 


= /i(fl® x - w - iqex /2 - ir ex / 2 )^,^ + 

This simultaneous equation set in Eq. (IA15I) is solved by the inverse matrix N(w) = [M(w)] _1 as 

(VM) = E JV ^'( W ) J dzV ^( z ) (Eo{z,u)). 

The incident light is assumed to be a monochromatic continuous wave laser as 

(E 0 (z,u])} = -E’o(z) <5(w - win). 


(A15) 


(A16) 


(A17) 


(A18) 
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Then, we rewrite Eq. (IA17I) as 

(Mw)) = j3 m (cj) S(u> - win), (Al9a) 

= ^2 N n, >*'(“) J d zV*i z ) E o{z). (A19b) 

A*' 

By performing Fourier transform, we can obtain an analytical expression of (6 Al (t)) under steady-state conditions as 

<M*)>=&(w (A20) 

By using this result, the simultaneous equations Eq. (IA12I) and Eq. (IA13I) can be analytically solved, and the population 
{w {t)bfj,(t)} and the correlation (6* are obtained under steady-state conditions. 


Appendix B: ENERGY CONSERVATION LAW 


In this section, we verify the energy conservation law, i.e., show that the incident light energy matches the total 
energy in the system. The total energy consists of the elastically and inelastically scattered photon energies, non- 
radiative decay energy, and dephasing energy. To calculate each type of energy, we follow the technique based 
on the input-output theory [3l|. First, we define the spectrum of the incident monotonic continuous wave laser 
Ei n (z,t) =£in S kl '*- UI °Q as 


S'in(w) = Re 


dr e iWT {£ in e ifeinZ }*£’ in = 4 S(u - w in ) 


(Bl) 


The spectrum of coherent components of the output radiation field on the transmission side is expressed as 


S c T M=Re 



e -+ - i 

dr -(£4(2o,0)) {Et(zo,t)) , 

7r J 


(B2) 


and that of one of the incoherent components, which corresponds to the photoluminescence (PL) spectrum (Eq.(4) in 
the main text), is expressed as 


S£(w) = Re 



e'“ T - + 

dr - (AE^(z o ,0)AE t (z o ,t)) 

7r 


(B3) 


Subsequently, we calculate the spectra of non-radiative decay and dephasing. In this work, we have employed the 
Markov approximation for the operators |_i = bj^b^+i and B M = bj^b^. Then, on the basis of the input-output 
theory, we can represent the output field as 

-^dampfW = — i Vlex G/j-l 
^phaseW = — W^ex B ^(t). 

In addition, within the Markov approximation, the number flux of the photons can be written as 

(B6) 

where the power of the light Q(z,t) in quantum optics is j28] 

Q(z,t) = 2£ 0 c{E\z,t)E(z,t )). (B7) 


(B4) 

(B5) 


Therefore, under steady-state conditions, the spectra of the non-radiative decay and dephasing are, respectively, 















TABLE I: Energy components in percentage for various film thicknesses and incident energies. 


L (nm) 

&j in (eV) 

Pc (%) 

Pic (%) 

-Pdamp 

(%) 

Pphase (%) 

-Ptotal (%) 

280 (nm) 

3.1860 

(eV) 

97.4253 

(%) 

2.4362 

(%) 

0.1499 

(%) 

-0.0114 (%) 

100.0000 (%) 

325 (nm) 

3.1860 

(eV) 

98.5195 

(%) 

1.4048 

(%) 

0.0808 

(%) 

-0.0052 {%) 

99.9999 (%) 

325 (nm) 

3.1936 

(eV) 

88.7898 

(%) 

10.5592 

(%) 

0.6716 

(%) 

-0.0207 (%) 

99.9999 (%) 

325 (nm) 

3.1981 

(eV) 

69.3738 

(%) 

28.6470 

(%) 

1.9680 

(%) 

0.0108 (%) 

99.9996 (%) 

325 (nm) 

3.1997 

(eV) 

58.2861 

(%) 

38.4866 

(%) 

3.2309 

(%) 

-0.0036 (%) 

100.0000 (%) 

325 (nm) 

3.2030 

(eV) 

94.8785 

(%) 

4.7041 

(%) 

0.4134 

(%) 

0.0039 (%) 

99.9999 (%) 


written as 

i r r°° p iWT 

SdanrpM = ^ ~ ( t )) 

A 4 

nj r r°° pi WT 

= -J- Re [5Z J o dT — (<W-i (0)<3 - m -i jA1 (r))j, 

A 4 

1 r r 00 pi^T 

<W>M = dr — ({^ hase (0)}t^ hase (r)) 

A 4 

r r . . r 00 p i<vr 

= ^ Re [^/ 0 dr ~ 

A 4 


(B8a) 


(B8b) 


where A = 2s 0 c/(fruj in ). 

In order to verify the energy conservation law, we have to establish the correctness of the following equation: 

/ oo poo poo poo 

dw hco S c (lo) + d U> hjJ S’ic(w) + / dw fuo Sdamp(w) + d LO hio S'phase(w) = filJi n £? n . (B9) 

-OO J— OO j — OO J —OO 

In TablejT] we show each component in percentage with respect to i.e., Pi = dw huj Si(co) / (hLo in £? n ) x 100. 

We see that the clephasing part P p hase sometimes becomes negative. Such a negative P p hase appears when the incident 
light energy is lower than the transverse exciton energy and the system can continuously absorb energy from the 
environment. This is discussed in detail in Ref. [3ll One can see that there is tiny mismatch between /total an d 100%. 
In order to check the validity of our numerical calculations, we examine the following quantity 

5 = X 100 (%) ‘ (B10) 

The quantity r] is the PL up-conversion efficiency from hui- m to the peak energy Hui 0 u t, where the output PL spectrum 
spreads between hui a and hutb- The quantity S roughly estimates the effect of the up-conversion on the total energy. 
For example, in case of the film thickness L = 325 nm and input energy hu> in = 3.186 eV, the up-conversion efficiency 
evaluated in the main text is ^ 73'212 Iv = 0.077 %, i.e., S = 0.0725 %. This is much larger than the maximum mismatch 
(~ 0.0004 %). Therefore, we can confirm that the mismatch is caused by numerical error and the energy conservation 
law holds true even though a portion of the incident light is up-converted. 


Appendix C: EFFECT OF THE NON-RADIATIVE DECAY BETWEEN EXCITON STATES 

In this section, we examine the effect of the non-radiative decay between exciton states on PL spectra. We perform 
calculations where non-radiative decay rate is zero and is comparable to dephasing rate hr ex = 0.2 meV. For instance, 
Fig.|3]shows the PL spectra in case of the input energy Hu> ln = 3.1981 eV and film thickness L = 325 nm. The red solid 
line represents the spectrum for fry ex = 0 meV, the green dashed line represents the spectrum for 7i7 ex = 0.01 meV 
(Fig. 2(a) in the main text), and the blue dotted line represents the spectrum for hj ex = 0.1 meV. In case of a larger 
non-radiative decay, the spectral shape of the green dashed line in Fig. [4] is slightly different compared with those of 
the other lines. However, the luminescent peaks appear at almost the same energies, even beyond the LT splitting. 
Therefore, we confirm that the mechanism of PL up-conversion does not change with non-radiative decay between 
exciton states. 
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FIG. 4: The PL intensity is plotted as a function of the output photon energy. The incident energy fk jj n in is tuned to the 
resonance energy of the exciton-light coupled mode m = 7 (3.1981 eV) in case of the film thickness L = 325 nm. The red 
solid line represents the spectrum in case of zero non-radiative decay rate Ir/ex = 0 meV. The green dashed line represents 
the spectrum in case of the non-radiative decay rate fr/ex = 0.01 meV in the main text in Fig. 2(a). The blue dotted line 
represents the spectrum in case of a larger non-radiative decay rate /r/ex = 0.1 meV. The numerical results indicate that finite 
non-radiative decay does not change the essential profile of the proposed up-converted PL effect. 
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